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The objectives of our investigation, namely "to characterize the atmospheric 
and directional effects on surface reflectance and vegetation index using the First 
ISLCSP Field Experiment (FIFE) dataset , develop new algorithms to obtain better 
AVHRR indices , and define possible improvements for future satellite missions", 
have been addressed in three separate, yet complementary studies. 

First, we have shown, from theoretical calculations, that visible and near- 
infrared reflectances combined linearly at optimum (one or two) viewing angles 
relate linearly to the fraction of photosynthetically available radiation absorbed by 
plants, fpar, can be used independently of the type of foliage and substrate, eliminate 
the effects of sub-pixel spatial heterogeneity, and improve the accuracy of the fp ar 
estimates when compared to the Normalized Difference Vegetation Index, NDVI. 

Second, we have demonstrated that NDVI, even though it is not a linear 
combination of radiances or reflectances, can be spatially integrated without 
significant loss of information from scales of 300 to 1000m. 

Third, we have successfully modeled AVHRR visible and near-infrared 
reflectances over the FIFE site, separating temporal and bidirectional components 
and determining the model parameters through an original iterative scheme. It 
appears that NDVI generated from the top-of-atmosphere reflectances normalized 
by the bidirectional effects (as determined in the scheme) is a better vegetation index 
than maximum NDVI. 

Details about the three studies are given below. 

1. Photosynthetically Available Radiation Abs orbed bv Plants 

We have completed a study of optimum combinations of visible and near- 
infrared reflectances for estimating the fraction of photosynthetically available 
radiation absorbed by plants, fpar- The results, published in the Proceedings of the 
5th International Colloquium on Physical Measurements and Signatures in Remote 
Sensing (Podaire et al., 1991; see appendix 1), strongly suggest that linear 
combinations of visible and near-infrared reflectances at specific viewing angles may 
be more accurate in predicting f par than vegetation indices that are non-linear 
functions of radiance or reflectance such as NDVI and the simple ratio, SR. For a 
sun at 60° from zenith in July at 45°N, for instance, the fp ar residual error is reduced 
from 0.058 to 0.033 (a factor of about 2). When using NDVI the minimum residual 
error is obtained for a nadir viewing, but when using linear combinations it is 
preferable to view the canopy at a 45° zenith angle. 

2. Spatial Integ ration of NDVI 

We have investigated the effect of spatial heterogeneity on the spatial 
integration of NDVI, which has concretized in an article published in IEEE 
Geoscience and Remote Sensing (Aman et al., 1991; see Appendix 2). For the 
considered sites, located in tropical West Africa and temperate France, and the scales 
analysed, 300 to 1000m, we have found that a strong correlation exists between 
NDVI calculated from average reflectances, MndvI/ an d NDVI integrated from 


individual NDVIs, Indvi- The relationship is almost perfectly linear, with a slope 
depending slightly on the variability of the vegetation cover. Effecting the scale 
change using Mndvi instead of Indvi does not introduce significant errors, 
especially when these errors are compared to those resulting from uncertainties in 
the relationships between NDVI and vegetation parameters, which are typically one 
order of magnitude higher. It is not excluded, however, that the principal scales of 
variability are smaller than the size of the pixels used in the calculations (20 to 30m), 
and other sites should be examined, as well as the variability of the vegetation below 
30 m, to conclude definitely and generally about the adequacy of using Mndvi i n 
biophysical parameterizations obtained at the local scale. 

3. Surface refle ctance and NDVI modeling 

Regarding our analysis of AVHRR data during FIFE, we have successfully 
modeled the channel 1 (visible) and channel 2 (near-infrared) reflectances. First, the 
data were carefully cloud-screened (Figs la, lb, lc. Id). This was done by examining 
the relationship between spatial standard deviation and average value over the FIFE 
site in channels 1 and 4 (thermal infrared) (Figs, la and lb), average values in 
channels 1 and 4 (Fig. lc), and spatial standard deviations in channels 1 and 4 (Fig. 
Id). Identifying on the plots the clear and cloudy cases as observed by an all-sky 

camera, we found that thresholds of 4 Wm'^sr^jim and 0.2 mWcm^sr'cm for the 
standard deviations in channels 1 and 4, respectively, provided the best cloud 
screening. All the data corresponding to channel 1 or channel 4 standard deviations 
above the threshold values were eliminated. Figs. 2, 3, and 4 illustrate the quality of 
the cloud-screening. Second, the raw data was transformed into reflectance (Fig. 5), 
from which was generated NDVI (Fig. 6). In the transformation, we took into 
account the drift in time of the AVHRR sensor according to Whitlock et al. (1990). It 
is interesting to note that the large variations in reflectance over a few days (Fig. 5), 
which camouflage any long-term trend, especially in the near-infrared (Fig. 5b), are 
substantially reduced when the visible and near-infrared reflectances are combined 
into NDVI (Fig. 6). The seasonal cycle of the vegetation becomes apparent, with high 
NDVI values in summer and low values in winter. Third, we expressed top-of- 
atmosphere reflectance as the product of a constant ,Ro/ and normalized temporal 
and bidirectional functions, F and G, respectively. Fig. 7 gives the model equations. 
For G, we used the simplest parameterization possible: G is expressed as the sum of a 
Lambertian contribution, a term corresponding to isotropic scattering in a semi- 
transparent medium (vegetation canopy), and a term taking into account the effect 
of shadowing structures. For F, we used a Fourier series limited to the 5th order 
(higher harmonics did not improve the accuracy of the modeling). Fig. 8 gives a 
shematic description of the iterative scheme developed to determine the model 
parameters. The bidirectional and temporal components of the model are estimated 
successively by linear, multivariate regression, and after a number of iterations, 
typically 10, convergence is obtained (see Fig. 9). Applying the model to the FIFE 
data, we explain 95% of the variance in both the visible and near-infrared channels 
(Fig. 10). Table 1 displays the model parameters, and Table 2 summarizes the 
comparison statistics between modeled and observed reflectances. The correlation 
coefficient is close to 0.98 for both channels and the standard deviation before 


regression is reduced from 0.45 to 0.10 and 0.21 to 0.01 in the visible and near- 
infrared channels, respectively. Fig. 11 shows RoF(t) as a function of time, and Fig. 12 
the resulting NDVI, NDVIt- Although RoF(t) and NDVI t are still contaminated by 
atmospheric effects (water vapor absorption and aerosol and molecular scattering), 
there are independent of solar and viewing geometries, which makes them useful 
to monitor the time changes of the earth's surface. In Fig. 11, for instance, a 
minimum of reflectance in the visible and a maximum in the near-infrared are 
observed during summer , and these features were not apparent in the raw data (Fig. 
5). In Fig. 12, the effect of vegetation stress around Julian day 220 (resulting from a 
lack of rain) is well depicted, but this stress could not be detected in the raw NDVI 
time series (Fig. 6). A strong correlation exists between observed NDVI and modeled 
NDVIt (Fig. 13), but the standard deviation remains rather large (12.5%) of the 
average NDVI value. Figs. 14 and 15 compare maximum NDVIs over 11 day periods 
(without and after cloud-filtering) and NDVIt, respectively. The correlation 
coefficients are large (above 0.98), but significant biases exist, especially when 
unfiltered data are used. It appears, for the dataset analysed, that the simple 
procedure of averaging cloud-filtered NDVI data yields better results than 
maximum NDVI (Fig. 16 and Table 3). 
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Upscale Integration of Normalized Difference Vegetation Index: The 

Problem of Spatial Heterogeneity 
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Table 1. Top-of-atmosphere reflectance model parameters. 
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Table 3 . Comparison statistics between maximum NDVI (without cloud-filtering) 
and modeled NDVIt ( case 1)/ between maximum NDVI (after cloud-filtering) 
and modeled NDVIt (case 2), and between average NDVI (after cloud-filtering) 
and modeled NDVIt (case 3 ). 
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Standard Deviation, Channel 1 (W/m2 sr Jim) 


30 



Average Radiance, Channel 1 (W/m2 sr (urn) 


Fig. la. Standard deviation of NOAA-9 AVHRR channel 1 (visible) radiance over 
the FIFE site versus average radiance. 
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Fig. lb. Standard deviation of NOAA-9 AVHRR channel 4 (thermal infrared) over 
the FIFE site versus average radiance. 
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Fig. lc. Average NOAA-9 AVHRR channel 4 (thermal infrared) radiance over the 
FIFE site versus average channel 1 (visible) radiance. 
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Fig. Id. Standard deviation of NOAA-9 AVHRR channel 4 (thermal infrared) 
radiance over the FIFE site versus standard deviation of channel 1 (visible) 
radiance. 
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Fig. 4. Time series of average NOAA-9 AVHRR channel 4 (thermal infrared) 
radiance over the FIFE site during 1987. (a) raw data; (b) cloud-filtered data. 
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Fig. 5. Time series of cloud-filtered average NOAA-9 AVHRR reflectances over the 
FIFE site during 1987. (a) channel 1 (visible); (b) channel 2 (near-infrared). 
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Fig. 7. Top-of-atmosphere reflectance model. Reflectance is expressed as the 
product of a constant,Ro, and normalized temporal and bidirectional functions, 
F and G. In the equations, t is time, T is the length of the dataset, P is the number 
of observations, N is the number of harmonics for F, and indices i and j refer to 
harmonics and observations, respectively. 



Fig. 8. Iterative algorithm to determine the reflectance model parameters. 
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Fig. 9. Residual error between NOAA-9 AVHRR reflectance model and observations 
as a function of algorithm iterations . (a) channel 1 (visible); (b) channel 2 (near- 
infrared). The various curves correspond to the numbers of harmonics used to 
describe F, 4 to 7. 
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Fig. 10. Observed versus modeled NOAA-9 AVHRR reflectances. (Top) channel 1 
(visible); (Bottom) channel 2 (near-infrared). Five harmonics are used to 
describe F. 
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Fig. 11. Time series of modeled RoF(t), observed R(t), and observed R(t) normalized 

by G(0 s ,0v/ ^s-^v) over the FIFE site during 1987. Satellite sensor is NOAA-9 
AVHRR. (a) channel 1 (visible); (b) channel 2 (near-infrared). 
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Fig. 12. Time series of Normalized Difference Vegetation Index calculated from 
modeled RoF(t), NDVIt(t), observed NDVKt), and observed NDVI(t) normalized 
by G(0s,0v, 4>s-<|)v) over the FIFE site during 1987. Satellite sensor isNOAA-9 

AVHRR. 




Fig. 13. Observed Normalized Difference Vegetation Index, NDVI, versus 
Normalized Difference Vegetation Index calculated from modeled R 0 F(t), 

NDVI t . 
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Fig. 14. Maximum Normalized Difference Vegetation Index over 11-day periods 
versus Normalized Difference Vegetation Index calculated from modeled RoF(t). 





Fig. 15. Maximum Normalized Difference Vegetation Index over 11-day periods 
(after cloud-filtering) versus Normalized Difference Vegetation Index calculated 
from modeled RoF(t). 
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Fig. 16. Average Normalized Difference Vegetation Index over 11-day periods (after 
cloud-filtering) versus Normalized Difference Vegetation Index calculated from 
modeled RoF(t). 
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ABSTRACT 

A useful parameter to estimate terrestrial primary productivity, 
that can be sensed from space, is the daily-averaged fraction of 
photosynthetically available radiation (PAR) absorbed by plants. 
To evaluate this parameter, investigators have relied on the fact 
that the relative amount of radiation reflected by a vegetated 
surface in the visible and near-infrared depends on the fraction of 
the surface covered by the vegetation and, therefore, correlates 
with absorbed PAR. They have used vegetation indices, namely 
normalized difference and simple ratio, to derive absorbed PAR, 
even from coarse spatial resolution sensors such as the Advanced 
Very High Resolution Radiometer (AVHRR) aboard the NOAA 
satellites. The problem with normalized difference and simple 
ratio is first, they are non-linear functions of radiance or reflectance 
and, therefore, cannot be readily applied to heterogenous targets, 
second, they are used in generally non-linear relationships, which 
make time-integrals of the indices not proportional to primary 
productivity, and third, the relationships depend strongly on the 
type of canopy and background. To remove these limitations, we 
propose linear combinations of visible and near-infrared 
reflectances at optimum (one or two) viewing zenith angles. 

Keywords: radiation, plants, primary production. 

1. INTRODUCTION 

Land primary productivity, or the rate at which materials from 
the atmosphere and soils are accumulated into biomass through 
photosynthesis, is of great importance. The reasons are numerous 
and have been discussed extensively in the literature (see, for 
instance, Ref. 1). Basically, the major portion of human food is 
provided by plants growing over land. Land primary productivity 
also affects the environmental context in which man and societies 
develop. In addition to its key role in sustaining human 
populations and structuring communities, land primary 
productivity governs to a large extent the seasonal oscillations of 
atmospheric carbon dioxide and impacts the water and energy 
available to the atmosphere. Anthropogenic changes, such as 
those linked to the destruction of major vegetation systems, have 
potential implications on climate. If we are to truly understand the 
interactions between the terrestrial biosphere and the atmosphere 
and their effects on climate, we need to know the geographic 
distribution and temporal variability of land primary productivity 
over the globe. 

To achieve this goal, satellite observations are essential. A 
promising technique for sensing primary productivity from space, at 
least in the case of light-limited situations, incorporates the fact 
that the growth rate of many plants is close to proportional to the 
rate at which radiant solar energy is absorbed by the foliage (Ref. 
2 ): 

PP =ef PAR PAR <d 

where PAR is the incident photosynthetically active radiation, 
practically the solar radiation reaching the canopy in the 0.4- 
0.7pm spectral interval, f PAR is the fraction of PAR intercepted by 


the canopy, and e is an efficiency factor for carbon fixation that 
depends slightly on plant type, temperature, and available soil 
water. For various canopies (mostly crops), t has been found to lie 
between 1.1 and 1.4 g C per M J of PAR (Ref. 3). 

Photosynthetically active radiation represents a nearly constant 
fraction of total insolation (e.g. Ref. 4), and total insolation can be 
retrieved accurately from satellite observations (e.g.. Refs. 5, 6, 7, 
8). Direct satellite estimates of PAR can also be obtained, as recent 
studies demonstrate (e.g.. Ref. 9). 

The absorbed fraction of PAR can be estimated from vegetation 
indices, the most commonly used being simple ratio, SR and 
normalized difference, ND. These indices are defined by: 



ND = 


R N + Ry 


( 2 ) 

(3) 


where Ry and Rtf are upwelling radiances in the visible and near- 
infrared (for instance radiances in channels 1 and 2 of the Advanced 
Very High Resolution Radiometer aboard NOAA satellites), 
respectively. Instead of radiances, reflectances are also used. That 
SR and ND are sensitive to f PAR results from the characteristic 
spectrum of sunlight reflected by leaves, which is distinct from 
that of sunlight reflected by soils. Chlorophyll pigments absorb a 
large fraction of the light which reaches them in the visible, but 
not in the near-infrared where scattering by the chloroplasts is 
effective. This is not the case of soils, whose reflectance increases 
more linearly with wavelength in the visible and near-infrared. It 
follows that the relative amount of radiation reflected by a 
vegetated surface in the visible and near-infrared depends on the 
fraction of the surface covered by vegetation and, therefore, 
correlates with f PAR . 

Several theoretical studies have predicted how SR and ND relate 
to f PAR (e.g., Refs. 2, 10, 11, 12, 13). Kumar (Ref. 2) suggested a near- 
linear relationship between SR and f PAR . Asrar et al. (Ref. 10) and 
Sellers (Refs. 11, 12) showed that f PAR should vary non-Iinearly 
with SR, but almost linearly with ND. Choudhury (Ref. 13) found 
that the relationships between f PAR and vegetation indices are 
generally non-linear. Soil reflectance changes, in particular, 
appeared to significantly affect the linearity of the relationships. 

Experimental studies (e.g.. Refs. 2, 14, 15, 16, 17) have also 
provided disparate results and, therefore, did not resolve the 
apparent theoretical controversy. Kumar (Ref. 2), for instance, 
observed that SR is linearly related to f PAR for sugar beet, which 
supported their theoretical analysis, yet Steven et al. (Ref. 14) 
reported an exponential relationship. In short, the observations 
indicate that the relationships between f PAR and vegetation 
indices depend strongly on the type of canopy and underlying 
surface, as pointed out by Choudhury (Ref. 13). 
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That the Interdependence of f PAR and vegetation indices is linear or 
non-linear is an important issue because linearity implies that the 
time-integral of SR and ND should also be linearly related to 
primary production (e.g.. Refs. 11, 12, 18). Even though serious 
doubts exist about the linearity of the relationships, especially 
when considering various ecosystems and climates. Tucker et al. 
(Ref. 19) and Coward et al. (Refs. 20, 21) have reported agreement 
between ND time-integrals obtained from AVHRR data over the 
Senegalese Sahel and North and South American biomes, 
respectively/ and published production rates. 

Using SR or ND for estimating f PAR , has several limitations. First, 
as mentioned above, the relationships are generally non-linear, 
which makes SR and ND time-integrals not proportional to 
primary production. Second, SR and ND are non-linear functions of 
radiance. Since vegetation is highly heterogeneous spatially, sub- 
pixel variability is likely to introduce uncertainties in SR and ND, 
particularly when the sensor spatial resolution is coarse (case of 
AVHRR). For such sensors, which have the advantage of frequent 
global coverage, applying relationships established for 
homogeneous canopies is not satisfactory. Third, satellite-derived 
SR and ND may reduce to some extent the effects of sensor 
calibration uncertainties and atmospheric interference, yet they 
depend on atmospheric composition, in particular aerosol and water 
vapor amounts, and viewing geometry (e.g., Ref. 22). Even time 
series of maximum AVHRR vegetation indices over a several-day 
period, which correspond to minimum atmospheric contamination, 
remain relatively noisy, and it has not yet been possible to identify 
whether the noise is due to residual variations in the atmospheric 
contribution or to variable directional surface properties. We need 
to address these limitations if we are to remotely sense f PAK from 
space accurately . 

Our objective, therefore, and the purpose of this paper, is to define 
optimum combinations of visible and near-infrared reflectances 
that: a) relate linearly to f PAR ; b) can be used independently of the 
type of foliage and substrate; c) eliminate the effects of sub-pixel 
spatial heterogeneity; and d) improve the accuracy of f PAR 
estimates when compared to SR and ND. 

Z METHODOLOGY 

Instead of using radiance ratios, we express fp^p as a linear 
combination of visible and near-infrared radiances or, 
equivalently, reflectances. This procedure, when applied to a 
coarse resolution sensor such as AVHRR, should eliminate or, at 
least, substantially reduce sub-pixel variability effects. Linear 
combinations of reflectances, known as "greenness" transformations, 
have been used for many years to study vegetation parameters, in 
particular by Hatfield et al. Ref. 15) and Asrar et al. (Ref. 17). 
These authors found that greenness obtained by combining 
reflectances measured by a Barnes Modular Multispectral 
Radiometer (MMR) in two visible and two near-infrared bands is a 
much more linear predictor of f P ^ R than simple ratio and 
normalized difference. They did not favor greenness, however, 
because of the smaller sensitivity of this index to f PAK and the 
strong dependence of the relationship between f PAR and greenness 
upon solar zenith angle and canopy geometry. If known (e.g., from 
theoretical calculations), the dependence upon solar angle or, more 
generally, radiation geometry, should not be a problem because 
solar and viewing angles can be determined exactly. The problem is 
to eliminate the effects of variable canopy geometry and soil 
reflectance in the relationships. Given a sun position, this may be 
possible for specific viewing angles. 

Our approach, therefore, is to simulate for varied soil and canopy 
parameters, namely leaf optical properties, soil reflectance, leaf 
area index (LAI), and leaf inclination distribution function (LIDF), 
above-canopy visible and near-infrared reflectances as well as f PAR 
and daily averaged (weighted by incident radiation) f PAR , f PAR . 
The simulations are performed which the SAIL canopy reflectance 
model (Ref. 23). Various radiation geometries (solar and viewing 
zenith angles, relative azimuth angle) are considered, as well as 
direct and diffuse fractions of incident solar radiation. The 
absorbed fraction_of PAR, f PAR , is computed as a function of solar 
zenith angle and f PAR as a function of latitude and season. From the 
reflectances, simple ratio and normalized difference are derived. 

We focus on f PAR since this parameter rather than f PAR is required in 


primary productivity models. In addition, since the sensors 
potentially useful to monitor land primary productivity from space 
are, or will be carried by heliosynchronous satellites (AVHRR on 
the" NOAA series, the POLarization and Directionality of the 
Earth Reflectance instrument, POLDER, on ADEOS, and the 
MODerate resolution Imaging Spectrometer, MODIS, on EOS) and, 
therefore, offer the possibility of viewing the same target under 
one or several geometries during a several-day period, we attempt 
to estimate f PAR from single or multi-angle combinations of visible 
and near-infrared reflectances. Indeed, the multi-angle approach 
is only suitable when the characteristics of the surface target do not 
change significantly over the several-day period. 

Thus, we regress Jpar at each latitude and month during the year 
against simple ratio, normalized difference, and visible and near- 
infrared reflectances. One and two viewing geometries are 
considered for the combinations of reflectances. The regression 
statistics, namely correlation coefficient, regression coefficients, 
and residual error of estimate are analyzed to determine the solar 
and viewing angles that minimize the effects of variable LIDF and 
soil reflectance. The improvement in the predicting power of the 
linear combinations is also assessed. 


3. RESULTS 

To illustrate our theoretical approach and show the promise of 
linear combinations, Figs. 1, 2, and 3 display selected results 
obtained with the SAIL model. The calculations were performed 
for LAIs of 0.01, 0.05, 0.1, 0.2, 0.5, 1, 2, and 5, erectophile, 
spherophile, and planophile canopies, soil reflectances of 0.1, 0.2, 
0.3, and 0.4, and typical leaf optical properties (reflectance and 
transmittance of 0.1 and 0.001, respectively, in the visible, and 0.45 
and 0.4 in the near-infrared). The soil reflectance was assumed 
white spectrally, and the various LIDFs were considered conjointly 
in the regressions. A US 62 standard atmosphere (Ref. 24) 
containing continental aerosols (Ref. 25) was overlying the canopy. 
We see in Fig. 1 (bottom) that the influence of the background on 
the /Tar versus ND relationships is substantial, especially at 
moderate LAIs, but is reduced considerably when using linear 
combinations of reflectances (Fig. 2, bottom). In this case, the points 
corresponding to a same LAI are generally more aligned with the 
best fit line. The result is a drastic improvement in the f PAR residual 
error. For the solar and viewing geometries of Figs. 1 and 2 , bottom, 
the residual error is reduced from 0.058 to 0.033. When using ND 
the minimum residual error is obtained for a nadir viewing (Fig. 1, 
top), but when using linear combinations it is preferable to view the 
canopy at a 45° zenith angle (Fig. 2, top). The f PAR residual error 
can be further reduced by combining linearly visible and near- 
infrared reflectances at two viewing zenith angles (Fig. 3). Using 
reflectances at nadir and 60° from zenith, for instance, reduces the 
f PAR residual error to 0.026 (Fig. 3, bottom). Smaller residual errors 
can even be obtained when the second viewing zenith angle is as far 
as possible from nadir ( Fig. 3, top). 

The above results, however, are only valid for a sun at 60° of zenith 
in July and at 45° latitude. For a sun closer to zenith, the minimum 
fpAR residual error is encountered at higher viewing zenith angles 
when using ND, for instance at 45° when the sun is at 30° from 
zenith. In the case of uni-angle linear combinations, the minimum 
at 45° (Fig. 3, top) moves to 60° when the sun zenith angle 
decreases to 30°. The picture is more complex with multi-angle 
combinations because of the many angular possibilities. In general, 
for a particular sun configuration^several viewing zenith angle 
pairs provide similar good results residual error around 0.020). 
For a sun at 30° from zenith, for instance, viewing at nadir, 15°, or 
30° from zenith and at 75° from zenith gives residual errors ranging 
from 0.019 to 0.021. The regression coefficients, however, are quite 
sensitive to the viewing geometries selected. 

4. DISCUSSION 

The results presented in section 3, although encouraging, should be 
interpreted with caution. No hasty generalization can be made at 
this point. First, the SAIL model has often showed weaknesses 
when compared to measurements; it does not predict a hot spot and 
is only appropriate for agricultural plants that form a layer-type 
canopy. More accurate canopy reflectance models may be used, at 
least to provide a reference. Second, the background reflectance 
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Fig. 1 Top: f PAR residual error as a function of viewing geometry for a 
sun at 60° from zenith in July at 45°N. In the calculations , f^ R 
estimated from normalized difference, ND, is compared_to f PAR 
obtained with the SAIL model. Bottom : scatter plot of fp AR versus 
ND for a nadir viewing . Variable soil reflectance results in points 
aligned rather perpendicularly to the best fit line (dashed line), 
especially at moderate LAIs f PAR and ND values around OS), which 
indicates that the relationship between f PAR and ND is not only 
non-linear, but also strongly depends on the type of background. 
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Fig. 2 Same as Fig. 1, but linear combinations of visible and near- 
infrared reflectances (single viewing geometry), p v and p N , 
respectively. Compared to Fig.l, the fp AR residual error is 

substantially reduced when viewing around 45° from zenith. Points 
corresponding to a same LAI, but different soil reflectances, are now 
more aligned along the best fit line. 


may vary with wavelength in the visible and near-infrared, as is 
the case with most soils (e.g., Ref. 27) and leaf litter (Ref. 28). The 
canopy may also be composed of living as well as dead leaves or, 
more generally, leaves of different optical properties. 
Calculations, therefore, should include more realistic situations. 
Still, our study strongly suggests that linear combinations at 
specific viewing angles may be much more accurate in predicting 
fpAR ^ an * nc K ces that are non-linear functions of radiances or 
reflectances. 

The relationships established theoretically, the gain in f PAR 
residual error when using preferential viewing geometries, the 
ability of the linear combinations to reduce soil and LIDF 
dependence, etc., remain to be verified using in situ measurements. 
Unfortunately, few data sets exist that contain concomitant /pa * and 
reflectance measurements at various viewing angles. The results, 
therefore, may not be statistically significant. A dedicated 
experiment to establish and verify the validity of the SAIL-based 
data fits, therefore, appears necessary. 

One should further emphasize that surface reflectances observed 
from space are inherently subjected to instrument noise and are 


contaminated by the atmosphere. Consequently, it will be necessary 
in the comparisons of the various estimators' performance to include 
the effects of instrument noise and atmospheric interference, which 
act differentially on simple ratio, normalized difference, and 
linear combinations. This can be done by simulating the top of 
atmosphere reflectances corresponding to the surface reflectances, 
correct those reflectances for atmospheric effects assuming typical 
atmospheric characteristics, and translate the effects of 
uncertainties in these characteristics into above-canopy reflectance 
uncertainties. The procedure is then to introduce the above canopy 
reflectance uncertainties in the regression datasets, as well as 
typical instrument noise. 

Our investigation should be viewed in the context of future 
spacebome radiometers, in particular MODIS on EOS and POLDER 
atmosphere reflectances for atmospheric interference, are 
particularly adapted to monitor /par and, therefore, primary 
productivity from linear (uni- and multi-angle) combinations of 
reflectances. This should lead, during the EOS era, to a better 
characterization of terrestrial primary productivity on a global 
scale. 
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Fig. 3 Same as Fig. 2, but two viewing geometries. Compared to Fig. 
2, the f PAR residual error is further reduced when reflectances at 
nadir and 60° from zenith are linearly combined. 
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Upscale Integration of Normalized 
Difference Vegetation Index: The 
Problem of Spatial Heterogeneity 


Angora Aman, Heremino P. Randriamanantena, Alain Podaire, and Robert Frouin 


Abstract — Spatial integration or radiometric parameters that 
describe the Earth’s vegetation cover is an important issue 
importa nt when studying global scale land-atmosphere inter- 
actions. The Normalized Difference Vegetation Index (NDVI) 
is currently used to characterize the land surface in terms of 
vegetation cycles and primary production. The NDVI is sensitive 
to fractional vegetation cover, canopy density, leaf architecture, 
and leaf physical state. A potential problem with NDVI, however, 
is that it does not depend linearly on radiance or reflectance. 
Since vegetation can be highly heterogenous spatially, and since 
the relationships between NDVI and vegetation parameters are 
established locally, using NDVI derived from coarse resolution 
sensors such as NOAA Advanced Very High Resolution Ra- 
diometer (AVHRR) may not be appropriate. The purpose of 
this paper, therefore, is to analyze the correspondence between 
NDVI calculated from average reflectances, Mndvi, and NDVI 
integrated from individual NDVI’s, Indm by simulating AVHRR 
data from high spatial resolution SPOT 1 Haute Resolution 
Visible (HRV) radiometer and Landsat Thematic Mapper (TM) 
data. For the considered sites, located in tropical West Africa and 
temperate France, and the scales analyzed, 300-1000 m, a strong 
correlation exists between the two types of index. The relationship 
is almost perfectly linear, with a slope depending slightly on the 
variability of the vegetation cover. Effecting the scale change 
using Mndvi instead of Indm does not introduce significant errors, 
especially w'hen these errors are compared to those resulting from 
uncertainities in the relationships between NDVI and vegetation 
parameters, which are typically one order of magnitude higher. 
Other sites should be examined, however, and the variability 
within TM and HRV Radiometer pixels quantified in order to 
conclude definitely and generally about the adequacy of using 
Mndvi in parameterizations obtained at the local scale. 


I. Introduction 

Global monitoring of the earth’s vegetation cover is im- 
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portant to understand land-atmosphere interactions and their 
effects on climate. This has been identified as a major task of 
the International Geosphere-Biosphere Program (IGBP) and 
its space component, the Earth Observing System (EOS) (see 
[1], [2]). Changes in land use and, hence, vegetation cover, 
directly impact surface water and energy budgets through 
plant transpiration, surface albedo, emissivity, and roughness. 
They also affect primary production and, therefore, the global 
carbon cycle. Since atmospheric carbon dioxide (CO 2 ) con- 
centration is mainly regulated by primary production, at least 
on a seasonal or annual time scale, characterizing vegetation 
dynamics is crucial to assessing the fate of the current increase 
in atmospheric CO 2 . 

Satellites can provide a spatially comprehensive view of 
land vegetation cover. The time frequency of the data, how- 
ever, may not be adequate to monitor phenological changes. 
This is the case of data from high spatial resolution sensors, 
such as the LANDSAT Thematic Mapper (TM) and the 
SPOT Haute Resolution Visible (HRV) radiometer, which have 
repetitive cycles of 16 and 26 days, respectively. Studies 
of vegetation dynamics, therefore, have been carried out 
using higher time frequency data from sensors with coarser 
spatial resolution, namely, the Advanced Very High Resolution 
Radiometer (AVHRR) aboard the NOAA satellites. 

The spatial resolution of AVHRR images (about 1.1 km 
at nadir) has been shown to be sufficient for investigating 
vegetation changes on subcontinental to global scales (e.g., 
[3]). However, because vegetation is highly variable spatially, 
a system with at least 500 m resolution will offer major 
benefits when compared to 1.1 km resolution system ([3], [4]). 

The AVHRR provides data in visible, near-infrared, and 
thermal infrared regions of the electromagnetic spectrum. 
There is a strong absorption of incident solar radiation by 
chlorophyll pigments in the visible, while scattering by the 
chloroplast in the near-infrared region leads to high reflectance 
values. Consequently, the Simple Ratio (SR) or the Normal- 
ized Difference Difference Vegetation Index (NDVI), which 
both exploit the differential reflectance properties of plants in 
the visible and near-infrared, have been used largely to monitor 
land vegetation.([5], [6]) The SR and NDVI are defined by: 

SR - NIR/VIS (1) 

NDVI = (NIR - VIS)/(NIR + VIS) (2) 

where NIR and VIS are satellite radiances or reflectances 
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in the near-infrared and visible, respectively. The SR and 
NDVI depend on vegetation parameters such as leaf greenness, 
fraction of vegetation cover, soil type, soil moisture (e.g., [7]), 
and these parameters vary according to climatic conditions. 
Maps of vegetation physical properties derived from SR or 
NDVI may be used as inputs into local, medium, or global 
scale models. It is highly desirable, therefore, to map them 
at compatible spatial scales. Using AVHRR-derived SR and 
NDVI, however, may not be appropriate because these indexes 
are nonlinear functions of radiances or reflectances. Since 
vegetation can be highly heterogeneous spatially, subpixel 
variability is likely to introduce uncertainties in SR and NDVI. 

In this paper, we address the problem of spatial integration 
of NDVI. For a given area, linear regressions are performed 
between mean and integrated NDVTs calculated from high 
spatial resolution data (HRV radiometer and TM). Mean 
NDVI is obtained from visible and near-infrared reflectances 
averaged over the area, whereas integrated NDVI is obtained 
by averaging all NDVI’s within the area. The scales considered 
range from 300 to 1000 m; they correspond roughly to the 
spatial resolution of present and future sensors dedicated 
to or suitable for large-scale vegetation studies.We discuss 
the influence of spatial heterogeneity and structure of the 
reflectance fields on the slope of the linear regressions. The 
results are presented for the Sudano-sahelian zone as well as 
temperate agricultural sites in France. 


II. Background 

Spatial integration of vegetation indices is necessary to 
study net primary production and land— atmosphere interactions 
at the global scale. The physical processes characterizing 
vegetation are often difficult to describe and are generally 
parameterized at the local scale using field measurements. 
The problem is to extend these parameterizations to larger 
scales when taking vegetation into consideration within Global 
Circulation Models (GCM’s). 

Many studies, experimental as well as theoretical, have 
shown that primary production depends on the fraction of 
photosynthetically active radiation intercepted by plants, or 
IPAR. Experimental measurements, on the one hand, suggest 
that IPAR depends linearly on NDVI ([8], [9]). Theoretical 
considerations, on the other hand, indicate that for IPAR 
values beyond 0.6 the relationship is highly curvilinear (e.g., 
[10], [11]). Therefore, the hypothesis of a linear relationship 
between NDVI and IPAR only holds for specific vegetation 
situations (e.g., well-mixed area, partially covered by vegeta- 
tion). 

Another important parameter used to characterize vegetation 
structure or architecture is the Leaf Area Index (LAI) ([12], 
[13]). It is related to many physical processes such as radiation, 
energy, or mass exchanges. This parameter can be successfully 
predicted from NDVI, as shown by [8] for wheat canopies 
in Arizona. The dependence between NDVI and LAI can be 
approximated by a linear relationship for LAI values lower 
than 3. For LAI values greater than 3, the dependence of LAI 
on NDVI is no longer significant. 


Thus, estimating IPAR or LAI from NDVI may not only 
be influenced by the nonlinear dependence of NDVI upon 
reflectance, but also by the nonlinearity of the relationship 
between NDVI and IPAR or LAI. Consider a landscape of 
surface st, composed of n elements. Each element T of this 
landscape is characterized by its reflectance R iy its vegetation 
index NDVI;, its leaf area index LAI;, its intercepted photo- 
synthetically active radiation IPAR,, and its surface The 
total reflectance R t of the surface s t is: 

n 

/?, = E Ri.sjst (3) 

so that 

n 

V is t = ^vis ; - Si / S( 

1=1 

n 

NIRt = Y. NIR * ' (4) 

i=l 

and the average NDVI over s t is: 

NDVI, = (NIR* - VIS t )/(NIRt + VIS,). (5) 


For this landscape, the actual LAI t and IPAR f are 
follows: 

defined as 

n 

LAIt = Y, LAI< • s,/s t 
1=1 

(6) 

n 

IPAR, = ]TlPAR t ’ *«/*• 
1 = 1 

(7) 

If we assume a linear relationship between NDVI and IPAR 
or LAI, (3) and (4) can be written: 

n 

LAI t = a y^NDVI< • Si/s t +0 

i= 1 

(8) 

n 

IPAR t = a'Y NDVI, ■ s,/s t +/?' 
1 = 1 

(9) 


where a, /3, a', and j3 f are constants determined at the local 
scale. We see from (5) and (6) that average radiometric 
measurements over s t could be used to estimate LAI t and 
IPAR* only if NDVI t , the average NDVI over s t , can be 
expressed as: 

n 

NDVIf = ^NDVI; • Si /s t . (10) 

i=l 

To illustrate the implications of this assumption, a theoretically 
simple case is considered: a site where all pixels are either 
composed of bare soil or dense canopy. All bare soil pixels 
have the same visible and near-infrared reflectances: VIS = 
0.165 and NIR = 0.188; and for dense canopy pixels: VIS = 
0.22 and NIR = 0.61 (these values correspond to typical 
in situ measurements at Ouango-Fitini (Ivory Coast). Fig. 1 
compares the resulting NDVIf with £ t ? =1 *NDVI;-Sj/s t when 
the fractional canopy cover varies from 0 to 1. In general, 
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Fig. 1. Comparison between Ikdvi and Mndvi when fractional vegetation 
cover varies from 0 to 1. Soil and vegetation radiances over Ouango-Fitini 
are used in the calculations. 

NDVI; is smaller than NDVIj • Si/s t . . The largest 

difference, 0.02 or 7%, is obtained for NDVI, values around 
0.3. This difference is not negligible; it would introduce similar 
relative errors on IPAR, (a' is close to 1), but larger errors on 
LAI, (a is typically equal to 2). In the limits (fractional canopy 
cover equal to 0 or 1), as expected there is no difference 
between NDVI, and]T" =1 NDVI, • s t /.s t . 

As mentioned above, the relationships of NDVI versus 
IPAR or LAI are curvilinear for plant canopies (bare soils) at 
high (low) IPAR or LAI values. Investigating the contribution 
of each surface element to the integrated value of NDVI, 
therefore, requires a knowledge of the spatial distribution of 
vegetation. The problem is fairly complex, especially when 
the vegetation cover is heterogeneous. 

III. Material 

A. Study Areas 

Three study areas, characterized by varied vegetation types, 
have been selected. The first two areas are located in West 
Africa, along a 550-km transect, covered by TM and HRV 
radiometer (Fig. 2(a)). This transect exhibits a large range of 
vegetation types (from woodland to shrub and grass steppe) 
and climatic conditions: the annual rainfall amount varies from 
480 mm at Bidi (Burkina Faso) to 1280 mm at Ouango-Fitini 
(Ivory Coast). The annual rainfall distribution is unimodal and 
the vegetation of this zone is organized according to the rainfall 
gradient. The third area corresponds to a temperate agricultural 
site located in Beauce (France). This site is dominated by 
winter crops, but forests and urban zones are also present. 
The study areas are described in more details below. 

1) Ouango-Fitini (Ivory Coast, 9 o 36’N-4°0r\V): This site is 
located in the northern limit of the Comoe National Park. The 
climate is tropical with a dry season from November to April 
and a rainy season from May to October. The main vegetation 
types are [14]: dense-tree savanna and woodland, scattered 
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Fig. 2. Location of the test areas in West Africa (a) and in Beauce, France 
(b) (Kong et ai (1988). 

tree savanna, and lowland and slope grass savanna, which 
represents 50, 30, and 10% of the total surface, respectively. 
The vegetation cover presents many large homogeneous zones 
and its seasonal evolution is related to alternating dry and 
rainy seasons. 

2) Bidi (Burkina Faso, 13 Q 55’N-2°30'W): This second site 
of tropical vegetation is located northwest of Burkina Faso, in 
the Yatenga Province. The rainy season is extremely short; it 
starts in July and ends in September. The main feature of the 
Bidi landscape is the scarcity of the vegetation cover except 
in the lowlands, where the tree cover (20-30%) includes tall 
trees and the grass cover reaches 90% at the end of the rainy 
season, but disappears during the dry season. 
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TABLE I 

Characteristics of the Data Used in this Study 


Dales of Data Acquisition 


Name 

Sensor 

Dates 

Ouango-Fitini 



wa 

LANDSAT/TM 

November 9, 1986 

wb, w2 


December 27, 1986 

wc, w3 


March 1, 1987 

wd, w4 

11 

May 4, 1987 

Bidi 



ba, bl 

LANDSAT/TM 

March 1, 1987 

be, b4 


May 4, 1987 

bb, be, bd 

SPOT-l/HRV 

October 4, 1987 

Beauce 



ca, cb 

SPOT-l/HRV 

May 1986 

cc, cd 


June 1986 

ce, cf, eg 

- 0 - 

October 1986 

Visible and Near-Infrared Wavebands of the Two Sensor Systems 

Satellite System 

Band Number 

Waveband (nm) 

LANDSAT/TM 

3 

630-690 


4 

760-900 

SPOT-l/HRV 

2 

610-680 


3 

790-890 


3) Beauce ( France , 48°25’N-1°30’W): The region is mainly 
agricultural. The main crops are wheat, barley, corn, and green 
peas, and they represent 45, 10, 8, and 7% of the total surface, 
respectively. The period of interest is May-June 1986, which 
corresponds to the winter crop maturity (Fig. 2(b)). 

B. Satellite Data 

Multitemporal LANDSAT/TM and SPOT/HRV radiometer 
digital images are the primary sources of data used in this 
study. Information from the visible and near-infrared channels 
of the two sensor systems are used to compute NDVI. The 
spatial resolutions of the TM and HRV radiometer are 30 
m and 20 m, respectively. Two subareas of 480x512 pixels 
are selected from Ouango-Fitini and Bidi at each date. Five 
subareas of the same size are selected from Beauce sites. 
For a given site, each subarea is representative of a type of 
landscape. Table I gives date and location of the images. 

IV. Methods 

As the goal of the study is to investigate the spatial 
integration of NDVI, the original images are degraded to 
simulate spatial resolutions of 300 to 1000 m. Although there 
are differences in band center wavelength and width between 
current orbital instruments, there is no significant differences 
between the NDVI derived from visible and near-infrared 
channels of sensors such as TM and HRV radiometer [15]. 

High spatial resolution data are often used for simulation 
purposes and are an excellent tool in preparing future space 
missions [4]. For test sites, such as ours, a 20-30-m spatial 
resolution allows one to observe homogeneous objects [16]. 
Different methods are used to coarsen the spatial resolution 
from high spatial resolution data. A method based on averaging 
reduces the noise by smoothing the frequency distribution [4]. 

It also reduces the occurence of both high and low values. To 
avoid underestimation of the predicted values, methods based 
on the modulation transfer function (MTF) between initial and 
desired data are often used [4]. We shall assume that at a 


spatial resolution coarser than 300 m, the MTF is a rectangular 
window, so that the MTF characterizing the original images 
has no effect in the following study. 

In our case, since the spatial degradation from original 
to coarser data is important (factor of 10 in resolution), we 
spatially average the data. From the original visible and near- 
infrared images of spatial resolution p , we generate a spatially 
degraded NDVI image of resolution H = hxp by aggregating 
h pixels in columns and p pixels in lines. The value of h is 
choosen to be a common multiple of the original resolutions 
of TM and HRV radiometer. 

If Xij and pij denote the visible and the near-infrared 
reflectances of a pixel (i, j), respectively, we define: 

h h 

on 

i=l j—\ 

and 

h h 

( 12 ) 

t=l j—\ 

where n = /i 2 . The mean NDVI is then 

Mndvi = {y - x)/(y + x ) (13) 

and represents the NDVI measured at the resolution H — hxp. 

The integrated value of NDVI at the resolution H, I NDV i, is 
defined as follows: 

h k 

Indvi = l / n EE z 'i (14) 

1 = 1 j — 1 

with 

Z ij (l/ij **“ “b Xij), (15) 

On each spatially degraded image, a method of linear re- 
gression is applied to determine the relationship between 
Mndvi and Indvi* The coefficients of the linear regression 
are determined by least square minimization. The regression 
equation of M N dvi versus I NDV i is: 

Mndvi = a Indvi + 0 (16) 

where 0 and ot are the intercept and slope obtained by linear 
regression, respectively. These two coefficients are given by: 

® DVI /^I NDVI (17) 

and 

0 = Mndvi - « Indvi (18) 

where r is the linear correlation coefficient between Indvi and 
Mndvi, and aM NDV1 and <ti * dv1 are standard deviations 
around the average values Mndvi and Indvi, respectively. 
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Fig. 3. Standard deviation <t Mn versus Mndvi over Ouango-Fitini on 
November 9, 1986 (a), Ouango-Fitjni on March 1, 1987 (b), Btdi on October 
4, 1987 (c), and Beauce in June 1986 (d). 


V. Results And Discussion 

A. Qualitative Description of Mndvi Variability 

The radiance measured by a coarse spatial resolution sensor 
such as AVHRR originates from different vegetation commu- 
nities and soils whose characteristic sizes are generally lower 
than the resolution cell. In order to investigate the significance 
of Mndvi provided by AVHRR, and understand the distribution 
of the objects inside each pixel, we have plotted the mean 
value of Mndvi in 1-km areas as a function of corresponding 
g'Mndvi- Figs- 3(a), (b), (c), and (d) show the relationship 
between <t M sdv , and Mndvi- 

At Ouango-Fitini, in November (biomass maximum), the 
values of Mndvi at 1 km resolution are relatively high (be- 
tween 0.5 and 0.6) and the values of <tm NDV] are low (less 
than 0.06). In March (end of the dry season), the decrease in 
Mndvi corresponds to the predominance of bare soil [17]. The 
highest values of M N DViare reached over islets of forest. The 
<7 m ND vt versus M N dvi values are plotted on Figs. 3(a) and (b) 
for March and November, respectively. 

At Bidi, in May (dry season), the values of Mndvi 
andcTM ND vi are low, between 0.15 and 0.2 and less than 
0.04, respectively, and are explained by the smallness of the 
vegetation cover and the brightness of the soil. In October 
(end of the rainy season), the highest values of MNDVI are 
associated with the highest values of ct Mndvi (Fig. 3(c)). A 
correlation coefficient of about 0.7 between M N dvi and <7m ndvi 
is computed, which is typical of heterogeneous vegetation. 
The highest values of Mndvi are in agreement with [18]. In 
contrast to Ouango-Fitini, at the end of the rainy season, the 
vegetation is not uniform. 

In Beauce, the behavior is different. In June, there is a global 
decrease of <t Mndvi as function of M N dvi (Fig. 3(d)). The 
highest M N dvi ’s correspond to low ct Mndvi ’s because they 
are associated with homogeneous fields. Here, the correlation 
coefficient between Mndvi and <tm ndv1 is about 0.8. 




Mean NDVI 


(b) 

Fig. 4. Histogram of Mndvi for Ouango-Fitini on November 9, 1986 and 
March 1,1987 (a) and for Beauce in June 1986 (b). 


Three kinds of behavior can be noticed in Figs. 4(a) and 
(b), which show the frequency distribution of Mndvi- (0 
at Ouango-Fitini the landscape is largely dominated by a 
homogeneous vegetation cover at the end of the rainy season 
(Fig. 4(a)); (ii) at Bidi, there is a considerable influence of the 
bare soil, so that the landscape remains heterogeneous even 
at the end of the rainy season (Fig. 4(a)); (iii) in Beauce, 
the associated histogram exhibits a left dissymmetry in June 
associated with the predominance of mature crops (Fig. 4(b)). 


B. Variability of Mndvi with Regard to Spatial Resolution 

Figs. 5(a), (b), and (c) show the relationship between 
heterogeneity and spatial resolution at various dates. M N dvi 
is calculated (using the method described in section 4) at 300, 
480, 600, and 960-m spatial resolutions. For each resolution, 
o Mndvi calculated and plotted against resolution. The stan- 

dard deviations corresponding to the original spatial resolution, 
i.e., 20 or 30 m, are also reported on the figures. 

For Ouango-Fitini, at the end of the rainy season (Novem- 
ber, Fig. 5(a)), the graph is flat and the ex m.m DV i l° w 0 ess 
than 0.25), showing that variability depends weakly on spatial 
resolution. In December (Fig. 4(a)), the graph exhibits a small 
peak at approximately 480 m. This behavior could be related to 
the regrowth of the leaves on the islets of forest [18] and could 
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Fig. 5. Standard deviation 0\i NDV i versus spatial resolution over (a) 
Ouango-Fitini, (b) Bidi, and (c) Beauce. 


be a manifestation of the dephasing cycle between herbaceous 
and ligneous stratas. 

For Bidi, in May (Fig. 5(b)), the graph is flat and the 
<7m ND vi l° w 0 n general less than 0.02). Since the vegetation 
is dry, the adjacent values of NDVI are strongly correlated, 
which induces little or no spatial variation from one pixel to 
another at coarse resolutions. At the end of the rainy season 
(October, Fig. 5(b)), the graph exhibits a high value at the 
initial resolution and remains flat between 300 and 480 m. It 
is difficult to relate this scale to a specific type of vegetation 
cover [19], but small patches of vegetation may be present in 
the landscape at this date. 

For Beauce, in May (Fig. 5(c)), the large decrease in 
^Mndvi between 30 and 300 m is due to the diversity of 
this agricultural region and to the mixing of active crops 
and bare soils or to poorly covered plots. In October, bare 
soils and harvested crops are predominant, and the decrease 
in variability between 30 and 300 m is small. 

The main conclusion is that the loss of spatial variability 


is larger between the original resolution of 20 or 30 m and 
300 m than between 300 and 1000 m, and that the variability 
is not significantly affected by spatial resolution in the range 
300-1000 m. 

C. Spatial Integration of NDVI 

It has been shown in [20] that a linear relationship ex- 
ists between AVHRR reflectances and Multi-Spectral Scanner 
(MSS) reflectances degraded to 1 km resolution. The problem 
is different for NDVI, because the equality between Mndvi and 
Indvi is algebrically false. In this section we quantify the 
relationship between Mndvi and I ND vi and study the influence 
of spatial variability on the relationship. 

1) Relationship Between M^dvi and I^ovt' The linear cor- 
relation coefficient r between Mndvi and Indvi at 1-km 
resolution (Figs. 6(a), (b), (c) and Table II) is high (above 
0.99). Comparable results are observed for other spatial res- 
olutions, namely, 300, 480, and 600 m for Ouango-Fitini, 
Bidi, and Beauce test sites, whatever the date of acquisition 
(see Table III for 300-m resolution). Thus, the high linear 
correlation between the two indices does not depend on time, 
spatial resolution, and land cover type. 

The slope a, on the contrary, varies with time, spatial 
resolution, and land cover type; but globally, it remains close 
to 1. Its variation with spatial resolution, weaker at Ouango- 
Fitini than at Bidi, will be analyzed below. 

The intercept j3 does not exhibit coherent variations. Its 
maximum value is lower than 0.02 over all test sites (Tables 
II and III). Furthermore the values remain close to zero at 
all simulated spatial resolutions. Globally, positive values 
correspond to low a values (/? > 0 when a < 1), 

2) Influence of the Spatial Variability on a: As reported, the 
correlation coefficient r between Indvi and Mndvi and the 
intercept /? do not present significant variations for all the 
scenes considered. The slope a is the only parameter whose 
variations with spatial variability needs to be studied. For 
a given spatial resolution H , a linear relationship exists 
between a\i NDVI and <ti ndvi with a very significant correlation 
coefficient of 0.999. This relationship can be written: 

*m ndvi = a(H) <r lNDV1 + b(H). (19) 

Combining (11a) and (12) yields: 

a(H, t) = r[a(H) + b(H) /<7 Indv , ( H , 0] (20) 

where t represents time (cri NDVl may vary with time for a 
given test site). 

On the one hand, (13) indicates that if H is kept constant, 
then a will depend only on <t I ndvi . Eq. (13) predicts that a 
increases when cri NDVI decreases, and conversely. This agrees 
with the observed variability of I NDV i (Section IV). Moreover, 
the results are valid for all spatial resolutions considered (see 
Tables II and III for 1000 and 300-m resolutions). 

On the other hand, the variations of <xi NDVI and a with 
respect to H when t is fixed are weaker (Figs. 5(a), (b), and 
(c)). When the spatial resolution coarsens from 300 m to 1 km, 
the coefficient a{H) decreases globally from 1.004 to 0.973 
at Ouango-Fitini, and from 0.980 to 0.946 at Bidi, while it 
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TABLE II 

Mean of Indvi (Eflnd)), Its Standard Deviation (^ind), Mean of U mV i ( E ( Mn 6 )), Its Standard Deviation (<7 Mnd ), Root Mean Square Error 
(RMSE) Between I NDV i and Mndvi, Correlation Coefficient (r), Slope (a), and Intercept (J) at 1-km Resolution for All Test Sites and Dates 


Resolution: 1 km 
Sites £(Ind) 

<T[nd 

£(Mnd) 

^Mnd 

RMSE 

r 

a 

t 3 

Ouango-Fitini 








w2 

0.376 

0.040 

0.379 

0.040 

0.00328 

0.998 

1.009 

-0.00101 

w3 

0.324 

0.037 

0.327 

0.036 

0.00378 

0.998 

0.983 

0.00817 

w4 

0.595 

0.055 

0.596 

0.055 

0.00179 

1.000 

1.008 

-.00366 

wa 

0.559 

0.017 

0.560 

0.017 

0.00086 

0.999 

1.014 

-0.00749 

wb 

0.372 

0.037 

0.374 

0.036 

0.00244 

0.999 

0.993 

0.00439 

wc 

0.337 

0.034 

0.339 

0.033 

0.00284 

0.998 

0.995 

0.00335 

wd 

0.569 

0.036 

0.569 

0.035 

0.00092 

1.000 

3.009 

-0.00480 

Bidi 









bl 

0.205 

0.029 

0.201 

0.026 

0.00552 

0.996 

0.904 

0.01554 

b4 

0.179 

0.009 

0.178 

0.009 

0.00067 

0.999 

0.969 

0.00501 

ba 

0.201 

0.034 

0.197 

0.030 

0.00587 

0.997 

0.903 

0.01528 

bb 

0.482 

0.045 

0.478 

0.042 

0.00557 

0.997 

0.950 

0.01940 

be 

0.511 

0.048 

0.506 

0.046 

0.00639 

0.999 

0.956 

0.01679 

bd 

0.527 

0.034 

0.521 

0.033 

0.00599 

0.998 

0.976 

0.00716 

be 

0.176 

0.011 

0.176 

0.010 

0.00066 

0.999 

0.976 

0.00383 

Beauce 









ca 

0.350 

0.048 

0.351 

0.048 

0.00334 

0.998 

1.023 

-0.00721 

cb 

0.305 

0.051 

0.304 

0.052 

0.00218 

0.999 

1.015 

-0.00539 

cc 

0.622 

0.036 

0.623 

0.036 

0.00404 

0.994 

0.997 

0.00308 

cd 

0.579 

0.051 

0.579 

0.052 

0.00296 

0.998 

1.010 

-0.00536 

ce 

0.001 

0.028 

0.008 

0.025 

0.00762 

0.994 

0.916 

0.00676 

cf 

0.010 

0.017 

0.015 

0.016 

0.00490 

0.995 

0.904 

0.00531 

eg 

0.026 

0.017 

0.030 

0.016 

0.00469 

0.995 

0.904 

0.00661 


increases from 1.04 to 1.08 in the Beauce region. Examining 
Figs. 5(a), (b), and (c) in detail, more complex variations 
corresponding to those of the slope a can be observed. The 
fluctuations of b(H) are more irregular and do not present any 
significant trend. In general, absolute values of b(H) remain 
lower than 10“ 3 , except for Beauce, and b(H)/<ri NDVI is about 
10 -2 . The variations of both a(H) and b(H)/crj NDV] are of 
about 10 -2 , for all land cover types and states. In summary, 
these variations are weaker than the variations of a and ai NDVI 
with respect to time. 

3) a Simulation : In order to further test the significance of 
the slope a simulated using (13), this slope is compared with 
that obtained from the linear regression between M NDVI and 
Indvi (Tables IV and V). Simulated and observed values of a 
are also plotted versus time on Figs. 7(a), (b), and (c). 

On average, there is no significant difference between 
observed and simulated slopes, except for Bidi and Beauce 
in November 1986 and October 1986, respectively, where 
it exceeds 0.045. Elsewhere, the maximum errors between 
simulated and observed a values are lower than 0.025 for 
all the studied test sites (Tables IV and V). 


species is not weighted by the fraction of the surface it covers. 
The resulting error can be estimated by the root mean square 
difference (RMSE) between the M N dvi and Indvi- 


RMSE = 


n 

1 fn ^ (Mkdvi 

i-l 


Indvi) 3 


1/2 


( 21 ) 


where n corresponds to the number of pixels at the resolution 
considered, 

This RMSE was calculated at the different dates, spatial 
resolutions, and test sites studied in this paper (see Tables 
II and III). The results indicate That the maximum error 
introduced when assimilating I NDV j to Mndvi is less than 
6.4 10" 3 . 

In order to interpret the RMSE values obtained, a limited 
development analysis is carried out. Each NIR and VIS is 
expressed as: 


NIR, = NIR + nir, 

VIS< = VIS + vis, (22) 


C. Accuracy of the Spatial Integration where m and m arc , h£ m£an va|u£S for a giv£n spa(iaI 

When computing M NDV i, the contribution of each vegetation resolution. Using these expressions to relate I ND vi to M NDV i 
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TABLE III 

Mean of Indvi (£(Ind), Its Standard Deviation (<r Ind ), Mean of ^/NDviC^Mnd)), Its Standard Deviation (<7 M nd)> Root Mean SQUARE Error (RMSE) 
Between /ndvi and M ndvi, Correlation Coetticient (r), Slope (o), and Intercept (,3) at 300-m Resolution for All Test Sites and Dates 


Resolution: 300 m 
Sites £(Ind) 

^Ind 

^Mnd 

^Mnd 

RMSE 

r 

a 

,3 

Ouango-Fitini 








w2 

0.391 

0.047 

0.392 

0.047 

0.00197 

0.999 

1.000 

0.00116 

w3 

0.319 

0.052 

0.321 

0.052 

0.00260 

0.999 

0.997 

0.00278 

w4 

0.575 

0.062 

0.575 

0.062 

0.00151 

1.000 

1.006 

-0.00271 

wa 

0.562 

0.024 

0.562 

0.024 

0.00077 

1.000 

1.006 

-0.00311 

wb 

0.387 

0.041 

0.388 

0.041 

0.00162 

1.000 

0.995 

0.00279 

wc 

0.327 

0.041 

0.328 

0.041 

0.00263 

0.998 

0.997 

0.00210 

wd 

0.574 

0.045 

0.574 

0.045 

0.00118 

1.000 

1.007 

-0.00351 

Bidi 









bl 

0.203 

0.035 

0.200 

0.033 

0.00444 

0.998 

0.938 

0.00948 

b4 

0.178 

0.011 

0.178 

0.011 

0.00057 

1.000 

0.981 

0.00301 

ba 

0.207 

0.051 

0.203 

0.048 

0.00576 

0.998 

0.940 

0.00884 

bb 

0.483 

0.061 

0.480 

0.059 

0.00483 

0.999 

0.967 

0.01218 

be 

0.511 

0.067 

0.507 

0.065 

0.00500 

0.999 

0.979 

0.00679 

bd 

0.527 

0.053 

0.522 

0.052 

0.00512 

0.999 

0.981 

0.00540 

be 

0.177 

0.015 

0.177 

0.015 

0.00061 

1.000 

0.982 

0.00276 

Beauce 









ca 

0.350 

0.082 

0.350 

0.084 

0.00386 

0.999 

1.013 

-0.00427 

cb 

0.305 

0.082 

0.304 

0.082 

0.00260 

1.000 

1.006 

-0.00234 

cc 

0.621 

0.072 

0.622 

0.072 

0.00455 

0.998 

1.002 

-0.00031 

cd 

0.578 

0.074 

0.579 

0.074 

0.00292 

0.999 

1.005 

-0.00274 

ce 

0.000 

0.042 

0.004 

0.041 

0.00525 

0.997 

0.961 

0.00387 

cf 

0.010 

0.030 

0.012 

0.029 

0.00375 

0.997 

0.958 

0.00304 

eg 

0.025 

0.030 

0.028 

0.029 

0.00338 

0.997 

0.963 

0.00316 


leads to (23) at bottom of page. Limiting the development 
of the previous equation to the second order in reflectance 
yields: 

f ndvi = Mmdvi(1 + A + B -f e) (24) 

where 

A = (Var(NIR) - Var(VIS))/(Nffi? - VlS*) 

B = (Var(NIR) + Var(VIS) + 2 Cov(NIR, VIS)]/ 

(NlR + VIS) 2 (25) 


and e is the residu. The variance and covariance functions, 
Var and Cov, are given by: 

n 2 

Var(NIR) = 1/n^ nir, 

i=l 

n 2 

Var(VIS) = l/n£ visi 

t=i 

n 

Cov(NIR,VIS) = 1/n^nir, * vis{. (26) 

t=i 

The terms A and B in (14) depend in a complicated way 
on the mean, variance, and covariance of the visible and near- 
infrared reflectances characterizing the pixels within the coarse 


n 

Indvi = 1 /n ^2 [NlR - VIS + ni ri - vis,] / [MR + VIS + n\Ti + vis,] 

1=1 

n 

Indvi = Mndvi 1 /n ^ [l + (nir,- - visj)/(NIR - VIS)] / [l + (nir ; + viSi)/(NIR + VIS)]. 

i=i 


( 23 ) 



334 


IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 30, NO. 2, MARCH 1992 





Fig, 6. Mean NDVI (Mndvi) versus Integrated NDVI (Indvi) on May 4, 
1987 at Ouango-Fitini (a), on October 4, 1987 at Bidi (b), and on October 
4, 1987 in Beauce (c). 


spatial resolution areas. SinceNIR is generally greater than 
VIS, the sign of A is essentially linked to the difference 
between Var(NIR) and Var(VIS). The term B can also be 
positive or negative depending on the sign and relative im- 
portance of Cov(NIR, VIS). A priori , there is no reason for 
the various parameters to combine such that ^4 + 5 = 0. The 
natural variability of the land cover and the reflectance values 
encountered in the real world, however, may minimize A + 5. 





Fig. 7. Time evolution of simulated and observed values of the slope a 
for Ouango-Fitini (a), Bidi (b), and Beauce (c). 

Let us now examine (24) using actual data. For the three 
sites and all dates, the absolute value of e remains small, 
typically 5 10“ 4 (see Table VI). Therefore, (24) may be 
justifiably written: 

Indvi = Mndvi (I + A + B). (27) 

The mean absolute values of Mndvi ( A + B) are small 
(less than 6.15 10 -3 ) for all sites and dates, which confirms 
the results displayed in Figs. 6(a), (b), and (c) and Tables 
II and III. The small difference between Indvi and Mndvi 
may be explained by the values of the parameters in the 
expressions of A and B (see Tables VI and VII). In general, 
Var(VIS), Var(NIR), and Cov(NIR,VIS) are^ much less than 
NIR 2 — VIS 2 and, a fortiori, (NIR -F VIS) 2 . The maximum 
value of Mndvi {A 4- B ), however, is 6.9 10“ 2 in the case of 
Bidi and Ouango-Fitini and 1.7 in the case of Beauce. This 
is due to pixels containing bare soil and vegetation (Ouango- 
Fitini) or bare soil (Beauce in Octobe r). When there is no 
vegetation, the difference NIR 2 - VIS 2 (i.e., the denominator 
in the expression of A) is small, yielding a high A value. The 
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TABLE IV 

Observed (a obs ) and Simulated (a simu ) Slope, with the Corresponding 
Correlation Coefficient (r), Standard-Deviation and Intercept 
(fi) for All Sites and Dates at 1-km Spatial Resolution 


Sites 

** obs 

Q simu 

r 

Ol 

a 

Ouango-Fitini: 





w2 

1.009 

1.000 

0.998 

0.0396 

-0.00101 

w3 

0.983 

1.001 

0.997 

0.0373 

0.00817 

w4 

1.008 

0.993 

0.999 

0.0548 

-0.00366 

wa 

1.014 

1.037 

0.999 

0.0172 

-0.00749 

wb 

0.993 

1.003 

0.999 

0.0369 

0.00439 

wc 

0.995 

1.004 

0.997 

0.0336 

0.00335 

wd 

1.009 

1.005 

0.999 

0.0355 

-0.00480 

Bidi: 






bl 

0.90401 

0.94853 

0.995 

0.0289 

0.01554 

b4 

0.96922 

0.96618 

0.999 

0.0093 

0.00501 

ba 

0.90304 

0.94871 

0.996 

0.0336 

0.01528 

bb 

0.95041 

0.94875 

0.998 

0.0451 

0.01940 

be 

0.95622 

0.94874 

0.998 

0.0484 

0.01679 

bd 

0.97575 

0.95002 

0.998 

0.0344 

0.00716 

be 

0.97557 

0.96283 

0.999 

0.0112 

0.00383 

Beauce : 





ca 

1.02273 

1.01026 

0.998 

0.0475 

-0.00721 

cb 

1.01487 

1.01684 

0.999 

0.0514 

-0.00539 

cc 

0.99699 

0.98654 

0.994 

0.0367 

0.00308 

cd 

1.01035 

1.01590 

0.998 

0.0513 

-0.00536 

ce 

0.91554 

0.95744 

0.993 

0.0278 

0.00676 

cf 

0.90443 

0.88800 

0.994 

0.0174 

0.00531 

eg 

0.90441 

0.88760 

0.995 

0.0173 

0.00661 


above cases, however, are isolated and cannot be considered 
as representative of the entire set of images. 

When the spatial resolution is 1 km, the average RMSE 
value for all cases considered is 0.0036, which corresponds to 
less than 1% of the average integrated NDVI value over all test 
sites and dates (see Table II). For individual sites and dates, 
the figures are similar, except for Beauce in Oct. (RMSE value 
reaching 0.0076). The 1% value is much smaller than that 
expected theoretically. In Section II we have indicated that for 
Ouango-Fitini the difference between Indvi and Mndvi should 
be about 6 % when Indvi is around 0.3. A possible explanation 
is that distinction between bare soil and vegetation cannot be 
clearly made even at the local scale, or that the principal scales 
of variability are smaller than the size of the TM and HRV 
radiometer pixels. 

It is necessary to compare the 0.0036 value to the errors 
linked, on one hand, to the NDVI itself, namely those resulting 
from corrections of atmospheric and surface directional effects, 
and, on the other hand, to the empirical relationships between 
NDVI and vegetation parameters such as IPAR and LAI. 
Atmospheric effects can induce errors of typically 0.02 and 
0.06 on the absolute NDVI value due to atmospheric water 
vapor and aerosols [21], respectively. The complete correction 
of these perturbing effects, especially aerosol scattering, as 


TABLE V 

Observed (a obs ) and Simulated (a sjmu ) Slope, with the Corresponding 
Correlation Coefficient (r), Standard Deviation (07), and Intercept 
03) for All Sites and Dates at 300-m Spatial Resolution 


Sites 

Os obs 

O' simu 

r 

<*! 

0 

Ouango-Fitini: 





w2 

0.999 

1.001 

0.999 

0.0476 

0.00116 

w3 

0.996 

1.001 

0.999 

0.0523 

0.00278 

w4 

1.005 

1.001 

0.999 

0.0620 

-0.00271 

wa 

1.006 

0.999 

0.999 

0.0245 

-0.00311 

wb 

0.995 

1.000 

0.999 

0.0416 

0.00279 

wc 

0.996 

0.999 

0.998 

0.0419 

0.00210 

wd 

1.006 

1.001 

0.999 

0.0454 

-0.00351 

Bidi: 






bl 

0.938 

0.963 

0.997 

0.0359 

0.00948 

b4 

0.980 

0.952 

0.999 

0.0115 

0.00301 

ba 

0.940 

0.965 

0.997 

0.0518 

0.00884 

bb 

0.967 

0.967 

0.999 

0.0619 

0.01218 

be 

0.978 

0.967 

0.999 

0.0670 

0.00679 

bd 

0.981 

0.966 

0.999 

0.0536 

0.00540 

be 

0.982 

0.957 

0.999 

0.0157 

0.00276 

Beauce: 





ca 

1.013 

1.008 

0.999 

0.0828 

-0.00427 

cb 

1.006 

1.008 

0.999 

0.0820 

-0.00234 

cc 

1.001 

1.002 

0.998 

0.0723 

-0.00031 

cd 

1.005 

1.005 

0.999 

0.0741 

-0.00274 

ce 

0.961 

0.976 

0.997 

0.0425 

0.00387 

cf 

0.957 

0.951 

0.996 

0.0304 

0.00304 

eg 

0.963 

0.952 

0.997 

0.0308 

0.00316 


well as surface directional effects is not possible. The empirical 
relationships between NDVI and vegetation parameters also 
exhibit significant errors. The accuracies of the associated 
physical parameterizations (e.g., [9], [22]) are often an order 
of magnitude greater than the RMSE values reported in this 
paper. For example, [9] found a RMSE of 0.053 in the linear 
relationship between NDVI and IPAR they established from 
field measurements. Thus the errors introduced when using 
Mndvi instead of Indvi are one order of magnitude lower than 
the residual errors after radiometric corrections of NDVI and 
those induced by the physical parameterizations. We conclude 
that using NDVI derived from coarse spatial resolution sensor 
(e.g., NOAA/AVHRR) data in relationships determined at 
the local scale will not introduce significant errors in the 
estimated biophysical parameters, at least for the situations 
selected and until more accurate relationships between NDVI 
and biophysical parameters are established. 

VI. Conclusion 

This study, achieved on two West African tropical sites and 
one French temperate site, shows that, for a spatial resolution 
lower than 1000 m (300—1000 m), the integrated NDVI, Indvi 5 
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TABLE VI 

Mean Values of Various Terms of the Second-Order Limited Development of I N dvi versus Mndvi at 1-km Resolution (see (24) ) 


Site 

£(Pr*) 

£(Var(VIS)) 

£(Var(NIR)) 

E(Cov) 

E(V IS) 

£(NIR) 

E 0) 

Ouango-Fitini: 
w2 -0.00234 

0.00008 

0.00056 

0.00013 

0.071 

0.160 

0.0000129 

w3 

-0.00257 

0.00013 

0.00072 

0.00013 

0.105 

0.207 

-0.0001526 

w4 

-0.00087 

0.00007 

0.00041 

-0.00003 

0.063 

0.251 

-0.0000079 

wa 

-0.00054 

0.00003 

0.00027 

0.00001 

0.075 

0.268 

-0.000012 

wb 

-0.00186 

0.00008 

0.00051 

0.00011 

0.082 

0.181 

-0.0000903 

wc 

-0.00147 

0.00018 

0.00075 

0.00019 

0.111 

0.225 

-0.0003377 

wd 

-0.00034 

0.00005 

0.00029 

0.00001 

0.066 

0.243 

-0.0000481 

Bidi 

bl 

0.00378 

0.00165 

0.00093 

0.00098 

0.258 

0.388 

0.0003063 

b4 

0.00049 

0.00039 

0.00043 

0.00036 

0.276 

0.396 

0.0000086 

ba 

0.00392 

0.00175 

0.00108 

0.00106 

0.260 

0.387 

0.0002967 

bb 

0.00439 

0.00036 

0.00021 

0.00008 

0.085 

0.242 

0.0000814 

be 

0.00553 

0.00043 

0.00020 

0.00009 

0.083 

0.254 

-0.0000002 

bd 

0.00546 

0.00041 

0.00013 

0.00007 

0.080 

0.254 

0.0001004 

be 

0.00047 

0.00041 

0.00046 

0.00038 

0.282 

0.402 

0.0000070 

Beauce 

ca 

-0.00058 

0.00022 

0.00040 

-0.00008 

0.076 

0.159 

-0.0001627 

cb 

0.00103 

0.00021 

0.00024 

-0.00006 

0.084 

0.158 

-0.0001823 

cc 

-0.00076 

0.00019 

0.00054 

-0.00012 

0.049 

0.214 

-0.0004332 

cd 

-0.00052 

0.00015 

0.00046 

-0.00007 

0.057 

0.215 

-0.0001040 

ce 

-0.00615 

0.00004 

0.00013 

0.00007 

0.060 

0.061 

-0.0004636 

cf 

-0.00398 

0.00003 

0.00009 

0.00005 

0.061 

0.063 

-0.0002916 

eg 

-0.00392 

0.00005 

0.00012 

0.00007 

0.064 

0.068 

-0.0001610 


TABLE VII 

Values of Various Terms of the Second-Order Limited Development of Indvi versus M ND vi at 1-km Resolution. 
They Correspond to the Parcel Where the Difference Between Indvi and Mndvi is the Most Significant (see (24)) 

Site 

Pr* 

Var(VIS) 

Var(NIR) 

Cov 

E(V IS) 

£(NIR) 

£ 

Ouango-Fitini 
w2 -0.02978 

0.00006 

0.00102 

0.00004 

0.060 

0.142 

0.0015177 

w3 

-0.05810 

0.00009 

0.00215 

0.00026 

0.087 

0.167 

-0.0003752 

w4 

-0.01257 

0.00012 

0.00187 

-0.00039 

0.052 

0.277 

0.0004974 

wa 

-0.00716' 

0.00003 

0.00192 

0.00009 

0.075 

0.302 

-0.0000599 

wb 

-0.01694 

0.00006 

0.00096 

0.00004 

0.074 

0.179 

-0.0002013 

wc 

-0.03796 

0.00021 

0.00219 

0.00047 

0.100 

0.188 

-0.0012388 

wd 

-0.00560 

0.00002 

0.00081 

-0.00004 

0.046 

0.234 

-0.0003404 

Bidi 

bl 

0.07546 

0.00618 

0.00221 

0.00331 

0.239 

0.399 

0.0013984 

b4 

0.00829 

0.00072 

0.00034 

0.00047 

0.281 

0.411 

0.0001481 

ba 

0.06870 

0.00581 

0.00179 

0.00279 

0.241 

0.407 

0.0015708 

bb 

0.02688 

0.00115 

0.00018 

0.00011 

0.084 

0.270 

0.0003420 

be 

0.03234 

0.00118 

0.00052 

0.00055 

0.082 

0.271 

-0.0001524 

bd 

0.02198 

0.00069 

0.00027 

0.00032 

0.069 

0.255 

0.0008218 

be 

0.01153 

0.00097 

0.00065 

0.00077 

0.257 

0.386 

0.0000256 

Beauce 

ca 

-0.01882 

0.00026 

0.00118 

-0.00036 

0.073 

0.193 

-0.0000233 

cb 

-0.01430 

0.00026 

0.00080 

-0.00025 

0.079 

0.173 

0.0006012 

cc 

0.04867 

0.00128 

0.00181 

0.00062 

0.058 

0.205 

-0.0228064 

cd 

0.02356 

0.00070 

0.00020 

0.00006 

0.075 

0.219 

-0.0024748 

ce 

0.29578 

0.00005 

0.00019 

0.00010 

0.051 

0.045 

0.0020951 

cf 

0.52351 

0.00004 

0.00018 

0.00008 

0.061 

0.059 

0.0005623 

eg 

-1.72490 

0.00020 

0.00043 

0.00028 

0.063 

0.064 

0.0033479 


Pr=A/KDvi (A + B) 
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is linearly correlated to the spatial average of NDVI, M NDV i . 
The slope of Indvi versus M N dvi depends only on the structure 
state of the surface, which can be expressed by the spatial 
variability of the canopy cover, in terms of the standard 
deviation <t I ndvi . 

In general, the slope and intercept of the linear regressions 
are close to 1 and 0, respectively, and, hence, the error induced 
when Indvi is used instead of M N dvi is small (average RMSE 
of 0.0036). This error is often an order of magnitude lower than 
either the residual errors after radiometric corrections on NDVI 
or the errors induced by the current physical parametrizations. 
The small differences between I N dvi and M NDV i are shown 
to be directly linked to the small variance and covariance of 
the visible and near-infrared reflectances Var(VIS), Var(NIR), 
and Cov(NIR,VIS) when compared to the mean quantities 
NIR 2 — VIS 2 and (NIR + VIS) 2 . Theoretical considerations, 
however, suggest that much higher differences should be 
found. The discrepancy might originate from the fact that bare 
soil areas are not apparent, even at the local scale; but it is also 
possible that the spatial variability is concentrated at scales 
smaller than the size of the TM and HRV radiometer pixels. 

Nevertheless, we conclude that NDVI, even though it is 
not a linear combination of radiances or reflectances, can be 
spatially integrated without significant loss of information. 
Further validations, however, need to be carried out for other 
land cover types than tropical savanna or temperate crop to 
assess whether using Mndvt instead of Indvi is generally 
adequate. 

The use of remotely sensed data to estimate biophys- 
ical parameters of vegetation dynamic models may solve 
a large part of the problems linked to spatial integration 
and scale change. The question, however, is to find inverse 
relationships between remotely sensed data and the biophysical 
parameters. The NDVI, thanks to its ability to be linearly 
integrated and to relate to major canopy parameters, appears 
as a useful radiometric index for the global parameterization 
of vegetation-related processes such as primary production 
and heat and mass transfer between the land surface and the 
atmosphere. 
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